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Abstract 

We present the first measurement of the vortex free-energy order parameter at weak 
coupling for SU(2) in simulations employing multihistogram methods. The result shows 
that the excitation probability for a sufficiently thick vortex in the vacuum tends to 
unity. This is rigorously known to provide a necessary and sufficient condition for 
maintaining confinement at weak coupling in SU(N) gauge theories. 
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The vortex free energy (also known as magnetic-flux free energy) order parameter 
in gauge theories is defined as the ratio of the partition function in the presence of 
a topologically trapped vortex excitation (introduced by a singular gauge transforma- 
tion) to that without it. Its Fourier transform w.r.t. to the center (Z(N)) of the gauge 
group (SU(N)) defines the so-called electric-flux free energy which is rigorously known 
to provide an upper bound on the Wilson loop. These flux order parameters can char- 
acterize all possible phases of a (pure) gauge theory, and furthermore do this in terms 
of the behavior of the excitation expectation for a vortex. They were first considered 
in the study of gauge theories in though the use of the analogous quantities in 
statistical mechanics goes back much further 0. The idea that vortex configurations 
underlie confinement at weak coupling has a long history, and has been the subject of 
intense recent activity. (We refer to Ref. |§ for a review of recent developments and 
references to early and recent work.) 

In view of the physical significance of the magnetic-flux free energy, it may appear 
surprising that it has not been measured in simulations over the last twenty years. 
Accurate determination of (differences of) free energies in gauge theories, however, is 
well-known to be difficult. In fact, it is at first not quite clear how one should go 
about computing such totally nonlocal (lattice-length) quantities. We present here a 
computation for the group SU(2) based on multihistogram methods Qj. Such a method 
was recently used in Ref. to compute the free energy of a pair of Z(N) monopoles, 
a quantity related to the 't Hooft loop operator. Our result demonstrates that the 
excitation expectation for a sufficiently extended 'thick' vortex at large (3 is essentially 
unity. This is the feature responsible for maintaining the confining phase in SU(N) 
gauge theories even at weak coupling. 

We work on a rf-dimensional hypercubic lattice A of size Li x • • • x with periodic 
boundary conditions in all directions. We generally denote bonds by b, plaquettes 
by p, cubes by c, etc. The plaquette action is denoted by A P (U P ), where, as usual, 
Up = Ilbep Ubi the product of the bond variables around the plaquette; for the minimal 
(Wilson) action A P (U P ) = — /3RetrC/ p . The trace "tr" is defined to include a 1/N 
normalization. 

A coclosed set of plaquettes (2-cells) is a closed set of (d— 2)-cells on the dual lattice. 
Thus, in d = 3, it is a closed loop of dual bonds; in d = 4, a closed two-dimensional 
surface of dual plaquettes. For fixed //, u, let denote a coclosed set of plaquettes 
that winds through every 2-dim [/iz/]-plane of A, i.e. a topologically nontrivial plaquette 
set wrapped around the periodic lattice (d-torus) in the (d — 2) directions A ^ fi, v 
perpendicular to /i, v This is depicted in figure [TJ(a) , where the short lines represent 
the plaquettes in V, with the horizontal axis representing the x^, x v directions, and 
the vertical axis the remaining [d — 2) perpendicular directions. 

Define the partition function 
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where the plaquette action A P (U P ) is replaced by the 'twisted' action Apir^Up) for 
each plaquette of V^. Here the 'twist' t^ v G Z(N) is an element of the center. There 
are thus (N — 1) different nontrivial choices for t^ v . The trivial element t^ v = 1 is the 
ordinary partition function Z\(l) = Z^. 
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Figure 1: Stack of plaquettes carrying twist winding around the periodic lattice, (a) 
and (b) are equivalent sets. 

As indicated by the notation on the l.h.s. of ([!]), the exact position or shape of 
is irrelevant; the only dependence is on the presence of the Z(N) flux winding through 
each [/iz/]-plane. It is indeed easily seen that V^ u can be moved around and distorted 
by a shift of integration variables, but not removed; it is rendered topologically stable 
by winding completely around the lattice (figure 0(b)). By the same token introducing 
two twists, on V^ u and r' on V' in (|lj), is equivalent to introducing one twist 
t" = t^t' since V^ u and V' can be brought together by a shift of integration variables 
(figure 0). This expresses the mod iV conservation of the Z(N) flux introduced by the 
twist. Thus, for N = 2, any odd number of such (nontrivial) twists is equivalent to 
one, and any even number to none. 



Figure 2: Equivalent sets V reflecting mod N conservation of the twist (N = 2). 
The magnetic-flux free energy order parameter is now defined as 

exp(-F mg (v,)) = — 
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Generalizations of (|I])-(|2p may be considered by introducing sets V K \ for several or 
all of the \d{d — 1) possible distinct choices of planes [k\]. 

The twist amounts to a discontinuous (singular) SU(N) gauge transformation 
on the configurations in (HD with multivaluedness in Z(N) (so it is single- valued in 
SU(N)/Z(N))), i.e. the introduction of a in(SU(N)/Z(N)) = Z(N) vortex. The set 
V^u represents the topological obstruction to having singlevaluedness everywhere. ([!]) 
is then the partition sum for the system with a topologically stable vortex completely 
winding around the lattice; and (0) is the normalized expectation for the excitation of 
such a vortex. Hence, it is also referred to as the vortex free energy. 

Choosing, say, [jiu] = [12] in ([!]), we now drop the \iv subscript. One is interested 
in the behavior of (Q) in the large volume limit (in the van Hove sense), i.e. as the size 
of the lattice increases in any power law fashion, e.g. L M = 2 la >* for some fixed choice of 
positive exponents a M , integer I — > oo. Let A = LiL 2 be the area of each [12]-plane, and 
L = L 3 ■ • • Ld the lattice volume in the perpendicular directions. One is interested, in 
particular, in L > A. The twist introduces a cost in action localized on the plaquettes 
in V. This cost, proportional to L, may be lowered if there are configurations that 
contribute with finite measure in the integral @, and allow the flux introduced by the 
twist to spread in the two directions perpendicular to V, so that the action is closer to 
its minimum; in other words, if there is finite probability for exciting a 'thick' vortex. 

For sufficiently large lattices, there are then three possibilities (r^ 1): 

(a) exp(-F mg (r)) ~ exp(-a(/3, r) L) 

(b) exp(-F mg (r)) ~ exp(-/3c(r) ^) 

(c) exp(-F mg (r)) ~ exp(-cLe^ /3 ' r ) A ) 

In case (a) the magnetic flux stays focused in a thin vortex; this describes a Higgs 
phase. In (b) the flux can spread in a Coulomb-like fashion lowering the free-energy 
cost; this describes a massless Coulomb phase, where the long distance behavior is 
accurately given by weak coupling pertubative expansion. In (c) the gain in thickening 
the vortex is exponential; this characterizes the confinement phase. It is important to 
note that, in contrast to (a)-(b), only (c) gives a value which survives and in fact tends 
(exponentially) to unity for all ways of taking the thermodynamic limit as described 
above; this is the signature of the confinement phase. 

Since our computation below is for N = 2, we now write explicit formulae only for 
this case. The Fourier transform of @ w.r.t. Z(N) is known as the electric-flux free 
energy. For N = 2 this is simply: 

exp(-F c i) = T ex P( -i? m g (T)) = 1 - exp(-F mg (-l)). (3) 

r=l -1 
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Consider now a rectangular loop C in a [12]-plane. Then, for any reflection positive 
plaquette action, the Wilson loop obeys the bound 0: 



tr(C/[C]))< (exp(-F el )J A , (4) 

where Ac is the minimal area bounded by C. (f|) shows that confining behavior (c) for 
the vortex free energy implies area-law for the Wilson loop with string tension bounded 
from below by the excitation expectation for a vortex. So confining behavior for the 
vortex free energy is a sufficient condition for linear asymptotic quark confinement. 

Placing suitable constraints in the functional measure in ([!]) which forbid the spread- 
ing of flux across [12]-planes, thus eliminating the occurance of thick vortices, results 
in nonconfining behavior of type (a) above |7[] . In this case (|]) cannot tell us any- 
thing about the Wilson loop. To show loss of confining behavior for the Wilson loop 
itself in the presence of such constraints, one needs a lower bound on it which exhibits 
perimeter-law. This was recently proven for large (5 in ||. Thus the occurance of thick 
vortices is also a necessary condition for confinement at weak coupling. 

Our measurement of (0) for SU(2) was done by an application of a multihistogram 
method [ffl. From now on we restrict the form of the action to the Wilson action 

A P (U P ) = -(3tvU p , (5) 

which was used in the measurement. The basic quantity in our procedure is the density 
of states w(S) as a function of the total action 5* along the twisted plaquettes. This is 
defined as 

w(S) =~[[[dU b exp (p ]T tiU p ) 5(S + ]T trU p ) . (6) 

P (fy pev 

If w(S) is known, the partition function can be easily computed for any coupling /3y 
along V as 

Z(J3v)= fdSw(S)e-^ s . (7) 



In particular, we are interested in Z(/3), the untwisted, and Z(—j3), the twisted parti- 
tion function. The problem is that the dominant contribution for Z(j3y) comes from 
different regions of S, depending on /3y. Therefore one needs to know w(S) to a good 
accuracy in a wide range of S. A simulation done at a certain value of however will 
give accurate information on w(S) only in a narrow neighbourhood of {S)/3 V . The main 
idea of the Ferrenberg-Swendsen multihistogram method is to combine information on 
w(S) coming from simulations at different /3y's to obtain w(S) in a wide range of S 
accurately. This can be done by noting that for a given (3\> the probability distribution 
of S, P(S,f3v), goes as 

P (S , MK _A_ e -^, (8) 
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and that P(S, /3\>) can be directly measured by making a histogram of the action along 
V. In this way, any simulation at a certain fly gives an estimate for w(S), 



w(S)=P(S,(3 v ) e^ s Z(f3 v ). 



(9) 



These estimates coming from simulations with different /3y's (say Pi, ^v-Pk) can then 
be averaged with suitable (S dependent) weights to minimise the error in w(S) over a 
given range of S. This results in the following set of coupled equations: 



which can be solved by iteration starting from Z(f3 n ) = 1. To optimize the procedure, 
one needs a sufficient overlap between the P(S, j3 n ) distributions corresponding to suc- 
cessive /3 n 's. Since the distributions quickly become narrower with increasing lattice 
size, the number of simulations, K, also needs to be increased accordingly. This makes 
our measurement very expensive on large lattices. For the largest lattices we typically 
used K = 40 — 80 simulations with the /3 n 's equally spaced in the — (3 +(3 range. 

The result of the computation for (fj) is shown in figure [| We have performed the 
computation on lattices of equal linear size in all directions for three different values of 
(3. The lattice spacings are a = 0.165 fm, a = 0.119 fm and a = 0.085 fm for (3 = 2.3, 
(3 = 2.4, and (3 = 2.5, resp. 

Notice that, with the lattice size expressed in physical units, the measurements for 
different /3's fall on the same curve, as they should. This indicates that the universal 
curve has been reached, and will not change at larger beta. Also, the onset of the 
sharp rise around 0.7 fm is in the region of the finite temperature deconfining phase 
transition providing another indirect consistency check. 

The approach to unity for sufficiently large lattice size in figure ^ is striking. In 
comparison, for Coulomb-like massless behavior, an upper bound obtained by action 
minimizing within the spin-wave approximation gives ~ exp( — (3 (vr/2) 2 ) « 0.085 at 
(3 = 2.3. The points forming the upper part of the plot are well within the confinement 
region. The string tension values extracted from the vortex free energy in the confin- 
ing region are consistent with the values from heavy-quark potential calculations (see 
eg. ||), though still better precision in the measurement of the vortex free energy is 
required for precise quantitative comparisons. 

In conclusion, the result of our computation clearly demonstrates that the weighted 
expectation for the excitation of a sufficiently thick vortex in the vacuum tends to one. 
In this sense the vacuum can indeed be viewed as having a 'condensate' of thick long 
vortices. This is sufficient for maintaining confinement at large (3 in SU(N) gauge 
theories. As mentioned, rigorous results also show it to be necessary: were the behavior 
for (fj) exhibited in figure |3] not to occur, confinement at large beta would be lost. 






(10) 





6 



3 sj; 

x = 2.3 
o (8=2.4 
» /S = 2.5 

® 

o 



0.0 0.5 1.0 1.5 

lattice size (fm) 

Figure 3: SU(2) vortex free energy (Jfy vs. lattice size 
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